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INTRODUCTION 


i,»  »,'  >,»  I,' 


The  transport  and  deposition  of  entrained  particles  from  the  flow 
stream  onto  the  walls  which  contain  the  flow  have  been  under  investigation 
for  the  past  thirty  years.  Because  of  the  complicated  nature  of  the 
problem,  research  work  has  been  limited  to  rather  simplified  conditions. 

In  general,  the  transport  mechanism  is  strongly  dependent  on  the  particle 
size  and  the  characteristics  of  the  flow  field.  Among  the  early 
investigations,  Friedlander  and  Johnstone*  developed  a  theory  for 
relatively  large  particles  (>ly)  and  Lin  et  al.^  investigated  the 
molecular-size  (dissolved)  particles.  For  small  particles,  the  transport 
is  dominated  by  the  Brownian  diffusion  while  for  larger  particles,  the 
particle's  momentum  has  to  be  included  in  addition  to  molecular 
diffusion.  To  include  the  inertia  effect  of  the  larger  particles,  it  is 
suggested  that  they  only  need  to  diffuse  to  within  one  stopping  distance 
from  the  wall;  from  that  point  on,  the  particles  will  coast  to  the  wall  by 
virtue  of  their  momentum.  Since  the  stopping  distance  could  be  appreciable 
compared  to  their  sizes,  the  deposition  rates  are  far  in  excess  of  those 
resulting  from  pure  diffusion.  For  small  particles,  however,  the  stopping 
distance  is  negligibly  small  such  that  diffusion  dominates  all  the  way  to 
the  wall.  Beal-*  developed  a  model  that  combines  both  concepts  of  small  and 
large  particles  discussed  above.  His  method  predicts  the  deposition  of 
particles  ranging  from  molecular  size  to  10y  in  a  channel  or  pipe  for 
steady  turbulent  flow.  Menguturk  and  Sverdrup^  applied  Beal's-*  method  to 


* Friedlander ,  S.K.,  and  Johnstone,  H.F.,  "Deposition  of  Suspended 
Particles  from  Turbulent  Gas  Streams,"  Industrial  and  Engineering 
Chemistry,  49 ,  1957. 

2 

Lin,  C.S.,  Moulton,  R.W. ,  and  Putman,  G.L.,  "Mass  Transfer  Between 
Solid  Wall  and  Fluid  Streams,"  Industrial  and  Engineering  Chemistry, 
45,  1953. 
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Beal,  S.K.,  "Departion  of  Particles  in  Turbulent  Flow  on  Channel  or 
Pipe  Walls,"  Nuclear  Science  and  Engineering,  40,  1970. 

^Menguturk,  M.  and  Sverdrup,  E.F.,  "A  Theory  for  Fine  Particle 
Deposition  in  Two-Dimensional  Boundary  Layer  Flows  and  Application 
to  Gas  Turbines,"  J.  of  Engineering  for  Power,  104,  69-70,  1982. 
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the  prediction  of  fine  particle  deposition  in  steady  two-dimensional 
boundary  layer  with  application  to  the  gas  turbine  operation.  For  even 
larger  particles  (  >10 u) *  inertia  effects  dominate  and  the  approach  is  to 
derive  the  equation  of  motion  for  each  particle  and  follow  their 
trajectories  (Lagrangian  method).  Ganic  and  Rohsenow^  investigated  the 
large  liquid  drop  deposition  in  two-phase  pipe  flow.  The  equation  of 
motion  for  each  drop  was  developed  taking  into  account  drag  forces, 
buoyancy  forces,  gravity  forces,  lift  forces,  and  reaction  forces  due  to 
evaporation.  The  deposition  rates  predicted  by  this  Lagrangian  model  help 
explain  the  behavior  found  in  measurements  of  boiling  heat  transfer  with 
dispersed  flow. 

As  a  result  of  the  literature  survey,  it  was  found  that  no  published 
reports  deal  with  the  transient  aspect  of  particle  transport.  The  unsteady 
particle  transport  is  encountered  in  many  short  time-frame  processes.  The 
deposition  rates  of  particles  are  usually  much  higher  and  some  special 
characteristics  may  not  be  predicted  by  a  steady  state  analysis.  In  this 
report,  the  transport  and  deposition  of  particles  with  diameters  ranging 
between  molecular  size  and  10u  are  analyzed  for  flow  impulsively  started 
over  a  flat  plate.  The  flow  stream  entraining  the  fine  particles  is 
numerically  simulated  by  the  vortex  sheet  model  proposed  by  Chorin^.  The 
transient  particle  transport  equation  is  solved  by  the  Strongly  Implicit 
Procedure  developed  by  Stone^.  The  objective  of  this  study  is  to  obtain 
the  particle  deposition  rates  as  a  function  of  time,  particle  size  and  flow 


condi t i ons . 


^Ganic,  E.N.  and  Rohsenow,  W.M.,  "On  the  Mechanism  of  Liquid  Drop 
Deposition  in  Two-Phase  Dispersed  Flow,"  J.  of  Heat  Transfer,  101, 
228-294,  1979. 

^Chorin,  A..).,  "Vortex  Sheet  Approximation  of  Boundary  Layers," 

J.  of  Computational  Physics,  2  7 ,  428-442  ,  1978. 

^Stone,  H.L.,  "interative  Solution  of  Implicit  Approximations  of 
Multidimensional  partial  Differential  Equations,"  SIAM  J.  Numerical 
Analysis,  5,  530-558,  1968. 
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MATHEMATICAL  MODEL 


A.  Flow  Field  -  Impiusive  Flow  over  a  Flat  Plate 


A  vortex  sheet  method  proposed  by  Chorin"  is  adopted  for  the  flow 
simulation. 

The  governing  equation  is  the  unsteady  Navier-Stokes  equation  in  terms 
of  vorticity,  and  the  schematic  is  shown  in  Figure  1. 
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where  £  is  the  two-dimensional  vorticity  scaler,  p  is  the  kinematic 
viscosity  of  the  fluid.  (j  =  (u,v)  is  the  velocity  vector,  u  is  the 
velocity  component  in  the  x-direction  and  v  is  the  velocity  in  the  y- 
direction.  It  is  noted  that  equations  (l),  (2),  and  (3)  include  the 
boundary  layer  approximations. 

The  initial  and  boundary  condi tons  are 


u  =  u  .  v  =  0,  t  =  0 


u*u  ,  y+«  and  t>0 
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where  u,,,  is  the  steady  free  stream  velocity. 

The  governing  equation  and  the  initial  and  boundary  conditions  may  be 
non-dimens ionali zed  by  the  following,  with  L  being  the  length  of  the  plate, 

Z  =  -^p,  u  *  U/uoo,  v  =  V/ua>,  x  *  X/L,  y  =  y/L  and  t  =  CU“/L  (7) 


The  dimensionless  governing  equation  and  the  auxiliary  conditions  then 


become 


35-35-35  1  3"5L  _  “ 

—  +  u  —  +  v  —  =  - — ,  — —  Re,  =  - 

3t  3x  dy  ReL  3y2  L  U 
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u  =  1,  v  =  0,  t=0 


u+1 ,  y  ♦  00  and  t>0 


u  =  0 ,  v  =  0 ,  y  =  0  and  t>0 


The  bars  will  be  dropped  from  this  point  on. 

Only  one  parameter,  Re,  results  in  the  dimensionless  system.  It  is 
well  known  that  the  numerical  solution  becomes  more  difficult  as  the 
Reynolds  number  increases,  and  the  analytical  solution  is  intractable. 
Chorin0  indicated  that  the  mesh  width  6  in  a  finite  difference  scheme  must 
satisfy  the  following  condition  in  order  to  maintain  stability. 

<$Re  =  0  (1)  (14) 

Therefore  for  relatively  high  Reynolds  number  flow,  the  number  of  mesh 
points  required  to  obtain  a  solution  would  be  prohibitive.  Accordingly,  it 
' of  interest  to  develop  a  grid-free  numerical  method.  One  approach  which 
has  been  used  with  considerable  success  involves  simulating  the  fluid  field 
using  discrete  vortex  elements.  The  interaction  between  these  vortex 
elements  is  then  used  to  model  the  developing  flow  field.  This  method  was 
first  used  by  Chorin^  and  has  since  been  applied  by  numerous  researchers  to 
various  turbulent  shear  flows  with  surprising  success.  These  numerical 
simulations  have  depicted  the  existence  and  interaction  of  large  scale 
vortex  structures  in  qualitative  agreement  with  flow  visualizations. 
Quantitative  results  of  mean  velocity  and  turbulance  intensities  have  also 
been  in  reasonable  agreement  with  experiments.  These  successes  have  been 
especially  impressive  in  view  of  the  fact  that,  in  most  cases,  the 
simulation  has  been  restricted  to  a  two-dimensional  model. 

•  ft 

Chorin  has  demonstrated  that  the  vortex  method  overcomes  many 
difficulties  encountered  with  the  finite  difference  scheme  for  flows  in 
separated  regions  or  in  the  wake  of  an  object.  Especially  for  the 
recirculating  turbulent  flow,  the  constant  change  of  flow  direction  and 
velocity  gradient  prevent  the  grid-system  approach  from  obtaining  accurate 
solutions  due  to  accumulated  truncation  and  round-off  errors.  Th“ 

"discrete  vortex  method"  is  essentially  a  numerical  simulation  of  the 

^Chorin,  A.J.,  "Numerical  Study  of  Slightly  Viscous  Flows,"  J.  of 
Fluid  Mechanics,  57,  785-795,  1973. 
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process  of  vorticity  generation  and  dispersal;  therefore,  it  is  grid- 
free.  The  vorticity  in  the  fluid  is  grouped  into  vortex  "blobs"  and  their 

« 

motion  is  governed  by  both  the  random  displacement  using  computer-generated  ' 

pseudo-random  numbers  and  mutual  interaction  effects.  The  process  is  time- 
dependent  and  all  the  vortex  blobs  have  to  be  followed  at  all  time  while 
they  are  present  in  the  system. 


Thete  direct  numerical  vortex  methods  require  little  or  no  empirical 
information  for  computing  the  flow  field  and  more  importantly  they  give 
transient  instantaneous  flow  information.  This  gives  them  a  decided 
advantage  over  empirical  time  average  turbulence  models. 

£  Q 

The  vortex  sheet  model0  is  an  extension  of  the  vortex  blob  method.0 
It  overcomes  both  the  convergence  problem  near  boundaries  and  the  specific 
vortex  blob  structure  requirement  of  the  vortex  blob  method. 

The  solution  of  equation  (8)  is  composed  of  two  parts.  First  the 
convective  motion  of  the  vorticity,  i.e.,  solution  of  the  following 
equation, 


3£ 

3t 


+ 


0 


(15) 


The  above  equation  corresponds  to  the  inviscid  part  of  the  equation 
(8)  or  when  Re-*-00. 

In  the  absence  of  boundaries  by  partitioning  the  vorticity  £  into  a 
sum  of  vortex  sheets  of  vorticity  £^, 


£  = 


£• 


(16) 


Equation  (10)  can  be  integrated  to  give 


v(y)  =  /q  u(x’z)  dz 


1  3 


(17) 


and  equation  (9)  yields 


» 

i 

p 

i 


(18) 


u(x,y)  =  1 


r 

-  £(x,z)dz 


It  is  seen  that  if  £  (x,z)  is  known,  (17)  and  (18)  yield  a  and  v. 

Based  on  (16),  consider  a  collection  of  N  segments  S.  of  vortex 
sheets,  of  intensities  £. ,  i  =  1,***,N  .  These  sheets  are  segments  of 
straight  line  such  that  u  on  one  side  of  Sj  and  u  on  the  other  side  of  Sj 
differ  by  £..  S;  is  parallel  to  x  axis  of  length  h  and  center 
at  r.  =  (Xj,y.)  .  The  motion  of  each  vortex  sheet  is  determined  by  all 

other  sheets  in  the  flow.  If  (uj,  vj)  is  the  velocity  of  S; ,  then 
according  to  equation  (18) 


u . 

i 


-  T  *  jfyj 


(19) 


where 


d  . 
J 


l-( | x .-Xj | /h) 


and  thi; 


sum  E  is  over  all  Sj  such  that 


(20) 


y . >y .  and  x . -x  .  <h 
J  1  1  '  J  1 


(21  ) 


and  the  vertical  velocity  v;  can  he  approximated  from  (17)  by 


VL  *  ( I1“I2)/h 


ll  '  <*l  *  h/2>n  * 


l2  =*  (x.-h/2)y.  -I  S.dTy* 
i  i  J*  J  J  J 


where 


dt  »  1-1  x.  +  ^/2  -  x . | / h 
j  I  L  jl 


d  .  =*  l-|x  .-h/2-x  .  I  /h 
J  1  i  J1 


y j  =  < y L »y j ) 


The  sum  is  over  all  Sj  such  that  O^dt^l  and  the  sum  is  over  all 
Sj  such  that  0^d*^l  . 

Once  (ui,Vi)  is  known,  the  vortex  sheet  Sj,  located  at  (x?,y?J  at  time 
step  n,  will  move  to  (x?+^,y?+^)  at  time  step  n+1  by  the  first  order 
accurate  finite  difference  equation 


and  the  sum  E  is  over  all 
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where  At  is  the  time  step  size  and  (u?,v?)  is  the  velocity  of  S;  at 


time  step  n. 


The  second  part  of  the  solution  is  dealing  with  the  viscous  diffusion 


port  ion. 


h_  i  y  s 

3t  '  Re  3y2 


with  intial  condition  fj(0)  =  5(x,y,t=0)  . 


A  solution  of  this  equation  using  random  walks  can  be  obtained  as 
follows.  Distribute  over  the  x-y  plane  points  of  masses  and 
locations  r.  =  (x.,y.),  i=l,***,N  ,  N  large,  in  such  a  way  that  the  mass 
density  approximates  £(0)  .  Then  move  the  points  according  to  the 
following 


y .  +  n. 

J  i  i 


where  Hj  is  a  Gaussian  distributed  random  variable  with  zero  mean  and 
variance  2At/Re.  The  u=l  at  y>0°  and  v=0  at  y=0  are  automatically 
satisfied. 

The  solution  to  Equation  (8)  is  then  the  linear  combination  of  the  two 
partial  solutions  given  by  equations  (28),  (29),  (31),  and  (32). 


n  .  .  .  n 

x .  +■  (  At) u . 


n  /  .  v  n 
y.  +  ( At)  v .  +  n; 


Equations  (33)  and  (34)  represent  the  solution  of  Equation  (8)  for  a 
system  with  no  solid  boundary.  But  u=0  must  be  satisfied  in  the  presence 
of  a  solid  surface.  According  to  Chorin^,  this  is  achieved  by  the 
vorticity  generation  at  the  solid  surface. 

(i  =  (u  ,v  )  is  the  flow  obtained  from  the  solution  procedure 
o  o  o 

described  above,  therefore,  it  does  not  satisfy  u=0  at  y=0.  At  the  wall, 
the  effect  of  viscosity  will  be  to  create  a  thin  boundary  layer  adjacent  to 
the  wall,  the  total  vorticity  in  the  layer  per  unit  length  of  the  wall  is 

boundary  layer 

/wall  Wy  =  uo  (y=,0)  (35) 

One  has  to  create  a  vortex  sheet  of  strength  uq  (y=0)  per  unit  length 
of  wall.  This  vortex  sheet  is  then  broken  up  to  elements  and  allowed  to 
participate  in  the  subsequent  motion  of  the  fluid  according  to  Equations 
(33)  and  (34). 

Thus  on  the  wall,  points  Q  ,•••,  are  designated  such  that 
distances  Q^C^,  ^2^3’ ***^m-l°m’  e9ua3  h  .  At  eac^  point  Q.  ,  the  tangential 
velocity  uq(q.)  's  evaluated  using  equation  (19).  Then  an  original  big 
vortex  sheet  of  strength  2uq  (Q.)  is  created  at  Q. .  The  random  wall  scheme 
requires  a  large  number  of  vortex  sheet  elements  to  converge  properly. 
Before  this  original  vortex  sheet  is  allowed  to  move  out,  it  is  broken  into 
a  number  of  vortex  sheets  of  equal  intensity  and  their  total  intensity  is 
equal  to  2uq  (Q.)  .  Then  these  group  of  vortex  sheets  are  allowed  to  move 
according  to  Equations  (33)  and  (34)  the  next  time  step. 

B.  Particle  Transport 

A  fluid  stream  containing  particles  flows  over  a  surface,  particles 
are  transported  and  deposited  by  various  mechanisms  depending  on  their 


basically  by  molecular  and  turbulent  diffusion  if  the  particle  size  is  less 
than  10  p.  But  simultaneously  the  momentum  of  larger  particles  should  be 
included  when  they  approach  the  solid  surface.  According  to  Fick's  Law, 
the  particle  flux  normal  to  the  flat  plate,  G,  is  given  by 


_  r.  3m 

G  =  -PIH— 
3y 


(36) 


where  p  is  the  mixture  density,  m  is  the  particle  mass  fraction  and  D  is 
the  binary  diffusion  coefficient,  and  it  may  be  written  as 


D=  (Dg  +  Dt) 


(37) 


Dg  is  the  Brownian  diffusion  coefficient  due  to  laminar  molecular 
diffusion.  It  is  given  by  Einstein  formula^. 


n  -  KfiT 
DB  3*pd 


(38) 


T  is  the  mixture  absolute  temperature,  dp  is  the  particle  diameter, 
p  is  the  fluid  dynamic  viscosity  and  Kg  is  Boltzmann  constant. 

is  the  turbulent  contribution  for  particle  diffusion.  Unless  the 

particle  is  extremely  small,  D^,  is  several  orders  of  magnitude  greater  than 

,0.  As  suggested  by  Liu  and  Ilori^,  D  is 
B  l 


^Einstein,  A.,  Theory  of  Brownian  Motion,  E.  P.  Dotton  Co. 

^Liu,  B.Y.H.,  and  Ilori,  T.A.,  "Aeorsol  Diposition  in  Turbulent  Pipe 
Flow,"  Environmental  Science  and  Technology,  _8_,  1  974  . 
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£j  is  the  fluid  momentum  eddy  diffusivity,  p^  is  the  particle  density, 

T  is  shear  stress  at  the  surface.  The  second  term  is  introduced  for 
w 

particle  size  correction. 

Next,  the  particle  transport  equation  in  a  boundary  layer  is 

described. 

Basic  assumptions  adopted  are  as  follows: 

1.  The  particle  mass  fraction  is  small  enough  so  that  the  fluid 
properties  are  unaffected  by  the  presence  of  the  particles. 

2.  Particles  are  small  such  that  molecular  and  turbulent  diffusions 
are  the  only  transport  mechanism  in  the  boundary  layer. 

3.  The  interior  of  the  particles  are  included  near  the  solid  surface 
using  the  free  flight  theory  of  Friedlander  and  Johnstone*. 

The  geometry  and  the  system  coordinates  are  shown  in  Figure  2.  The 

conservation  of  particles  is  governed  by 
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Figure  2.  Particle  Transport  Mechanism  at  the  Boundary 
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The  boundary  condition  at  the  wall  needs  some  special  attention  in 
order  to  include  the  particle  inertia  effects.  Because  the  presence  of  the 
wall,  the  particles  behave  differently  from  their  diffusion-convection 
transport  in  the  interior  of  the  boundary  layer.  According  to  the  theory 
of  Friedlander  and  Johnstone^,  particles  of  relative  larger  diameter  coming 
that  close  to  the  wall  need  only  a  small  impetus  to  reach  the  surface. 

This  impetus  is  the  particle  inertia.  In  more  detail,  particles  transport 
toward  the  wall  occurs  by  diffusion  until  they  reach  the  stopping  distance 
from  the  wall.  Particles  are  considered  to  complete  the  final  phase  of 
their  travel  by  way  of  a  free  flight  (inertia  coasting)  which  takes  place 
under  the  momentum  imparted  by  molecular  Brownian  motion  or  turbulent 
fluctuations.  By  assuming  Stokes  flow  and  taking  into  account  the  particle 
radius,  the  particle  stopping  distance  may  be  calculated  from  the  following 
equation 
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where  S  is  the  stopping  distance  and  vg  denotes  the  free  flight  initial 
velocity.  As  Beal^  points  out  vg  may  be  determined  by  summing  the  Brownian 
velocity  due  to  molecular  motion  and  root-mean-square  normal  fluctuation 
velocity  at  the  stopping  distance.  Therefore, 
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The  stopping  distance  S  and  initial  velocity  vg  are  solved  by  iteration 
from  Equations  (44)  and  (45). 

The  particle  flux  at  the  stopping  distance  (y=S)  is 
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For  simplicity,  the  particles  crossing  y=s  line  are  assumed  to  be 
transported  to  the  wall  at  an  average  velocity  v^  which  appoximately  is 
equal  to  vg/2.  The  flow  of  particles  that  reach  the  wall  and  are  captured 
is  then  given  by 


Tppms 


where  P  is  the  fraction  of  the  particles  which  stick  on  impact.  The 

average  velocity  is  halved  because  there  is  equal  probability  of  particles 

being  thrown  back  toward  the  main  stream.  m  is  the  mass  fraction  of 

particles  at  the  stopping  distance.  Assuming  local  equilibrim  which  was 
•> 

tested  by  Beal  ,  then 
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The  above  is  the  required  boundary  condition  at  y=S. 

Next,  we  will  non-d imens i ona l i ze  the  transport  equation. 

Let  u=a/uoi,  v=v/vco,  Om/m^,  x=x/L,  y=y/L,  t=tuoo/L  ,  equation  (AO)  and  the 
auxilliary  conditions  become 
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In  matrix  form,  Equation  (53)  becomes, 
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In  the  above,  the  subscripts  i,j  denote  the  grid  point  locations  and 
the  superscript  N  denotes  the  time  step  number. 

The  initial  and  boundary  conditions  become 
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where  n=l  corresponds  to  t=0+,i=l  to  x=0,  i=M  to  x= 1 ,  j  =  l  to  Y=S/L  and  j=N 
to  y+°°. 

The  governing  equations  in  matrix  form,  with  five  nonzero  diagonal 
elements  as  the  coefficient  matrix,  may  be  written  as  follows 
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In  order  to  accelerate  the  solution  procedure,  matrix  [M]  is  modified  as 
matrix  [M+N],  which  consists  of  seven  nonzero  diagonal  elements.  Equation 
(64)  then  becomes 


[M+N]  [  0n+1  ]m+1  =  [  Q°]  +  [N]  [0n+1r 


(66) 


in  which  m  is  the  number  of  iterations  for  determining  the  column 
matrix  [0]  at  the  (n+1)  th  time  step.  The  modifier  matrix  [N]  had  to 
fulfill  the  requirnient  that  the  coefficient  matrix  [M+N]  can  be  factored 
into  the  product  of  a  lower  matrix  [L]  and  an  upper  matrix  [V],  each 
consists  of  three  nonzero  diagonal  elements  in  the  lower  and  upper 
portions.  The  purpose  of  this  modification  is  to  determine  the  column 
matrix  [0°  ^]  by  having  the  difference  of  the  (m+l)th  and  the  (m)th 
iterations  within  a  prescribed  tolerance.  If  the  column  matrix  [AO]  at  the 
(m+l)th  interation  is  defined  as 
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then  Equation  (65)  can  be  written  as 


[M+N]  [A0n+1]m+l  =  [ ©n ]  —  [  M ]  I  0n+1  ]m 


If  one  uses  the  lower  and  upper  matrices  given  the  [L]  and  [U], 
Equation  (67)  can  be  expressed  as 


[L][U]l  A0n*1]m+1  =  [Qn )  -  (M][0n+1]m 


The  use  of  matrices  (Lj  and  [U]  accelerates  the  solution  procedure  by 
successive  eliminations.  Once  the  column  matrix  [A0n+^]  at  (m+1)  th 
iteration  is  equal  to  or  less  than  the  prescribed  tolerance,  the  value  of 
[0n+^l  at  the  (m)th  iteration  become  the  solution  of  the  new  mass  fraction 
of  the  particle  at  (n+l)th  time  step. 

A  computer  code,  named  PARTRAN,  has  been  developed  based  on  the  model 
presented  in  the  above  MATHEMATICAL  MODEL  section. 


VERIFICATION  OF  THE  COMPUTATIONAL  SCHEMES 


A.  Flow  Field  Calculations 

The  vortex  sheet  method  was  checked  with  the  well  known  Blasius 

1 2 

solution  of  the  steady  state  flow  over  a  flat  plate  .  Due  to  the  nature 
of  the  random  vortex  sheet  method,  the  instantaneous  results  are  expected 
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to  show  some  fluctuations.  Therefore,  when  comparing  with  the  steady 
laminar  solution  of  Blasius,  only  time-averaged  velocity  profiles  are 
appropriate  for  comparions.  The  test  case  uses  the  following  input 
pa  ram- t - rs 

=  lo\  At  =  0.2,  h  =  0.2 

The  instantaneous  drag  coefficient  and  the  averaged  drag  coefficient 
are  plotted  in  Figure  4.  The  averaging  is  carried  out  every  twenty  time 
steps. 

Figure  5  shows  the  location  of  vortex  sheets  at  t=6.4.  It  shows  the 
boundary  layer  development  and  the  distribution  of  vorticity. 

Figure  6  shows  the  instantaneous  velocity  profiles  at  five  different 
downstream  locations  at  t= 6.4. 

Figure  7  displays  the  comparison  of  the  three  averaged  velocity 
profiles  with  the  Blasius  velocity  solution. 

Based  on  the  comparisons  of  drag  coefficients  and  velocity  profiles, 
it  may  be  claimed  that  the  vortex  sheet  method  is  doing  an  adequate  job. 
The  drag  coefficient  actually  represents  the  bulk  effects.  The  time 
averaged  drag  coefficient  goes  from  0.429/Re^^  for  the  first  twenty  time 
steps,  then  to  0.642/Re  the  second  twenty  time  steps,  and  to  0.66/Re 
the  third  twenty  time  steps.  The  Blasius  steady  state  drag  coefficient  is 
0.664/Re^2.  This  tells  that  the  flow  will  reach  steady  state  in  about 
thirty  to  fourty  time  steps  (6  to  8  dimensionless  times).  The  steady  stat 
drag  coefficient  predicted  by  the  vortex  sheet  method  is  almost  equal  to 
the  exact  solution.  The  variance  for  the  drag  coefficient  is  around 
0.03/Re1/2. 

The  velocity  profile  comparisons  show  the  same  trend  as  the  drag 
coefficient.  The  deviations  are  all  less  than  1"'. 


vigure  5.  Typical  Vortex  Sheet  Positions  in  Boundary  Layer 
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B.  Particle  Transport  Computations 


B.l  Comparison  with  Exact  Solutions 


In  order  to  compare  with  the  exact  solution,  one  of  the  boundary 
conditions,  Equation  (52),  has  to  be  changed  to  the  following 


0,  y  =  0 


The  Tr  number  is  therefore  dropped  out  of  the  system.  Physically  it  means 
that  wall  is  a  perfect  sink  for  the  particle;  whenever  a  particle  reaches 
the  surface  it  will  be  removed. 

With  Equation  (69)  replacing  Equation  (52),  the  flow  and  particle 
transport  equations  with  auxiliary  conditions  (8)~( 12)  and  (49)~(52)  are 
reduced  to  the  following  for  steady  state  condition 
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( 70)  and  ( 71 )  become 


f  "  +  l/2ff"  =  0 
C(0)  =  0,  f  *  ( °°) -*-1 
f(0)  =  0 
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The  exact  solution  to  the  above  set  is  given  in  (12)  also, 


O’(0)  =  0.332  Sc1/3Re1/2 


96(0) 

9y  average  between  0  and  x 


=  0.664  /Re/x  Sc 


Two  cases  were  selected  for  comparison.  The  results  nr  :  shown  in 
Table  1.  Again  the  solutions  from  the  computer  code  PAKTRAN  were  averaged 
over  twenty  time  steps 
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TABLE  1 


30(0)  /  3y  by  Eq.(79)  30(O)/3y  by  PARTRAN 


589.57 

2099.75 


592.86 

2050.5 


Based  on  Table  1,  the  good  agreement  shows  that  PARTRAN  Is  working 
correctly. 

B.2  Steady  State  Comparison 

This  section  describes  the  comparison  of  the  steady  state  solution  of 
PARTRAN  with  the  solution  of  Blasius  flow  coupled  with  the  steady  transport 
equation. 

The  governing  equations  are 
f'"  ♦  l/2ff"«0 
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f(0)  -  0 

30  30  m  _1  320 

3x  3y  Pe  3y2 

9.1,  <82) 

30 
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The  computer  program  BLASIUS  was  written  to  solve  the  set  of  equations 
(81)  and  (82).  The  Blasius  flow  is  solved  by  the  shooting  technique 
involving  fourth-order  Runge-Kutta  method.  The  steady  transport  equation 
is  solved  using  the  same  Strongly  Implicit  Procedure. 

The  comparison  between  the  results  from  PARTRAN  and  those  of  BLASIUS 
are  shown  and  discussed  in  next  section  -  Results  and  Discussion. 


SOLUTION  PROCEDURE 


i 

i 

Because  of  the  highly  fluctuating  nature  of  the  vortex  sheet  flow 
j  calculations,  proper  time-averaging  of  the  velocity  is  necessary  before 

using  it  for  particle  transport  calculation.  In  this  analysis,  since  there 
i s  no  feedback  of  the  particle  on  the  flow  modeling,  the  complete  flow 
field  was  calculated  and  plotted  first.  By  examining  the  flow 
fluctuations,  piecewise  time  averaging  was  performed.  The  averaging  period 
is  corresponding  to  the  fluctuating  period  of  the  vortex  sheet  flow.  This 
smoothed  velocity  field  is  then  entered  in  the  particle  transport 
calculation  using  Equation  (54).  This  procedure  substantially  reduces  the 
numerical  errors  in  the  finite  difference  procedure  solving  particle 
transport . 


RESULTS  AND  DISCUSSION 

The  problem  we  are  investigating  corresponds  to  the  physical  situation 
of  a  suddenly  started  impulsive  flow  of  a  fluid  containing  particles.  This 
impulsive  flow  has  a  uniform  velocity  of  uo  before  it  is  over  the  flat 
plate.  The  mass  fraction  of  the  particle  in  the  flow  before  the  flat  plate 
is  m^.  The  kinematic  viscosity  of  the  fluid  is  u  and  density  p.  The 
length  of  the  plate  is  L.  The  particles  in  the  flow  are  assumed  spherical 
with  diameter  dp  and  density  Pp.  The  binary  mass  diffusion  coefficient 
between  the  particle  and  the  fluid  is  D. 

It  is  noted  that  even  though  the  geometry  chosen  here  is  a  simple  2- 
dimensional  cartesian  system,  the  results  and  trends  obtained  can  be 
applied  to  impulsive  flow  inside  a  relatively  large  cylindrical  tube.  The 
approximation  is  valid  as  long  as  the  boundary  layer  thickness  is  small 
compared  with  the  tube  diameter. 

The  transport  equations  and  their  initial  and  boundary  conditions  are 
all  represented  with  dimensionless  variables.  The  results  if  the 
e a  1  ill  a t i ons  are  also  presented  in  terms  of  the  dimensionless  groups  in  the 


After  examining  the  dimensionless  governing  equations  and  their 
auxiliary  conditions,  three  dimensionless  groups  appear  in  the  system. 

These  are 

Re  -  u  L/v  (83) 

oo 

Sc  -  v/D  (84) 

Tr  «  2D/vfPL  (85) 

Therefore  the  solutions  of  the  problem  may  be  represented  as  a 
function  of  the  above  dimensionless  groups. 

The  selection  of  the  ranges  of  these  dimensionless  groups  is  based  on 
practical  interests. 

The  Reynolds  number  varies  between  10^  and  10^.  10^  is  about  the 

lower  bound  for  the  boundary  layer  theory  to  be  valid  while  10^  is  in  the 
transition  regime.  The  vortex  sheet  method  has  been  tested  to  handle  flow 
up  to  106  by  Chorin6,  but  its  ability  to  deal  with  the  turbulent  Elow  is 
untested.  Theoretically  there  is  no  limit  to  Reynolds  number  for  the 
vortex  sheet  model.  It  should  be  pointed  out  that  the  flows  containing 
particles  or  pollutants  are  likely  to  have  earlier  transition  from  laminar 
to  turbulent  flows  and  the  transition  may  last  longer  than  pure  flow.  The 
instantaneous  flows  predicted  by  the  vortex  sheet  method  or  other  vortex 
blob  method  inherently  include  those  fluctuations  which  may  be  due  to  the 
presence  of  the  particles. 

The  Schmidt  number  also  has  a  wide  range.  For  laminar  diffusion  of 
particles  of  the  molecular  size,  e.g.,  liquid  in  another  liquid,  or  gas  in 
gas,  or  gas  in  liquid,  the  Sc  varies  from  0.74  (Oxygen  in  air)  to  1630 
(Glycerol  liquid  in  water).  When  the  diffusing  particle  size  gets  bigger 
so  does  the  Sc  number.  For  laminar  diffusion  of  particle  of  lu  in  air 
stream  using  Equation  (38),  the  Sc  number  is  around  10^.  The  Sc  number 


will  be  smaller  if  the  turbulent  contribution  is  included  [See  Equation 
(37)].  Also  the  Schmidt  number  is  around  0.9  for  turbulent  diffusion  of 
molecular  size  material.  Based  on  the  above  discussion,  the  Sc  number  is 
set  to  vary  between  0.7  and  10^. 

The  transport  number  Tr,  according  to  Equation  (52),  is  involved  with 
D,  v^,  P,  and  L.  Therefore  potentially  Tr  will  have  quite  a  large  range. 
For  a  ly  particle  in  air  and  laminar  flow,  the  stopping  distance  is  on  the 
order  of  10~^m,  the  v^  is  10-^  m/s,  and  Tr  is  around  10-^  with  P  assumed 
about  unity.  Therefore  Tr  would  get  larger  with  small  sticking 
coefficients.  Tr  is  increasing  with  decreasing  particle  size  and  it  is 
usually  smaller  in  turbulent  flow  than  in  laminar  flow  for  the  same  size 
particle.  To  cover  most  cases  of  practical  interest,  Tr  is  given  a  range 
of  l  to  10  ^ . 

In  the  numerical  calculation,  At  was  set  between  0.1  and  0.2,  Ax  was 
fixed  at  0.025.  Ay  was  determined  based  on  the  following  equation 

Ay  =  0.00005  x  1000  /  /Re~  (86) 

The  above  equation  is  designed  to  accomodate  the  boundary  layer  size  change 
at  different  Reynolds  numbers  because  the  total  number  of  nodal  points  is 
fixed.  There  are  A0  streamwise  nodes  and  160  vertical  nodes. 

Based  on  the  stopping  distance  theory,  the  vertical  ( y-d i r ec t i on) 
nodes  should  begin  at  the  stopping  distance  (y=S)  rather  than  at  the 
wet  1 1  .  After  evaluating  the  order  of  magnitude  of  the  stopping  distance,  it 
was  found  that  the  maximum  stopping  distance  is  on  the  order  of  10  ”m 
according  to  Equation  (4 A).  Since  the  Ay  is  on  the  order  of  10-/4m,  it  is 
reasonable  to  neglect  the  exi stance  of  the  stopping  distance  in  the  nodal 
arrangement . 

The  deposition  flux  of  particles  on  the  plate  surface  is  given  by  the 
following  dimensional  form 
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But  we  will  present  the  deposition  flux  in  dimensionless  form 
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The  dimensionless  flux  is  termed  Nusselt  number  or  Sherwood  number. 

In  our  study,  Nu( t)  is  defined  as  the  average  flux  at  each  time  step.  It 
is  calculated  by  taking  an  average  over  the  total  Nu  (t,x)  at  those  40 
nodal  points  on  the  wall. 

First  we  will  investigate  the  effects  of  different  Reynolds  numbers  on 
the  deposition  of  particles.  The  system  parameters  are  kept  constant 
except  the  Reynolds  number  Sc  is  set  at  10^  and  Tr  at  10  The  Reynolds 
number  varies  from  10  to  10  .  The  instantaneous  average  Nussett  number 
Nu(  t) ,  the  dimensionless  deposition  rates,  is  plotted  in  Figure  8.  As 
mentioned  previously,  the  random  nature  of  the  flow  field  induces  a  similar 
fluctuating  particle  transport  rate.  Because  all  flow  fields  were 
generated  with  the  same  random  number  pattern,  the  random  patterns  for  the 
deposition  rates  are  more  or  less  identical.  The  differences  among  these 
transport  rates  are  the  magnitudes.  As  expected,  higher  Reynolds  number 
generates  higher  deposition  rates.  From  the  statistical  point  of  view, 
more  meaningful  information  for  random  fluctuations  is  obtained  by  time¬ 
averaging  the  instanteneous  values.  Figure  9  shows  the  corresponding  time- 
averaged  Nussett  numbers  Nu  (t)  .  The  averaging  was  performed  for  every 

five  time  steps,  i.e.,  one  dimensionless  time.  The  curves  are  fairly 
smooth  and  reveal  some  trends.  After  very  steep  drops,  they  all  reach 
steady  state  around  t=5.6  dimensionless  time.  It  is  noted  that  in  real 
time,  they  may  not  reach  steady  state  at  the  same  time  because  t=tu  /L. 
After  some  analysis,  it  was  found  that  the  steady  state  values  are  well 
correlated  by  the  following  equation 
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(Nu  )  steady  state 
t  «a . _ _ _ 

/Re 

Another  interesting  information  to  examine  is  the  accumulated  deposition 
W(t),  as  a  function  of  time,  which  is  plotted  in  Figure  10.  The  curves  go 
up  quickly  in  the  beginning  and  then  taper  off  as  deposition  rates 
dec  rease . 

A  useful  formula  also  was  formed  for  the  accumulated  depositions, 
W(t,Re  )  W(t,Re0) 

- — -  =  - - —  (90) 

/Re^  /Re^ 

The  steady  state  comparions  between  the  results  of  PARTRAN  and  those  of 
BLASIUS  are  shown  in  Figure  9.  The  dashed  lines  represent  the  predictions 
by  BLASIUS. 

The  variation  of  the  Sc  on  the  particle  deposition  rate  is  examined  in 
Figure  11,  where  Re=l0^,  Tr=10-^,  and  Sc  varies  from  l  to  10^*.  The 
transport  rates  increase  with  increasing  Sc  for  all  the  transient  cases. 
Therefore,  the  same  trend  extends  from  the  steady  state  transport  to  the 
unsteady  transport.  In  the  steady  state  Equation  (79)  suggests  that  Nu  is 

.  1  /O  ... 

proportional  to  Sc  '  .  The  dashed  lines  indicate  the  steady  state  values 

of  the  depositon  rates,  Nu  predicted  by  the  program  BLASIUS.  It  is 

s  •  s  • 

seen  that  PARTRAN  is  doing  a  good  job  for  Sc  of  1,  10,  and  10^.  For 
intermediate  Sc  numbers,  the  steady  state  transport  rates  predicted  by 
PARTRAN  are  much  lower  than  those  from  BLASIUS.  It  is  believed  the  results 
of  BLASIUS  are  correct.  The  errors  of  PARTRAN'S  predictions  are  considered 
primarily  due  to  the  truncation  errors  of  the  finite  difference  method, 
ihe  statistical  errors  of  the  flow  fields  predicated  by  the  vortex  sheet 
method  amplify  the  errors  of  the  finite  difference  scheme.  The 
fluctuations  shown  in  Figure  8  explain  the  errors.  The  predictions  of 
the  Nu(t)  before  t<  1.5  are  reasonably  good,  but  the  large  peak  at  t  =  2 
introduces  the  main  truncation  errors.  It  is 
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interesting  to  note  that  if  the  curve  were  extrapolated  after  t=2.5  with 
the  same  slope  as  that  before  the  slope  change,  the  prediction  would  match 
very  well  with  the  steady  state  solution  predicted  by  BLASIUS  as  shown  by 
the  chain  line. 

Figure  12  represents  the  effects  of  the  Tr  number  on  the  time-averaged 
transport  rates.  As  shown  in  the  figure,  the  smaller  the  Tr  number,  the 
larger  the  deposition  rates.  Since  Tr  is  the  ratio  of  boundary  resistance 
to  interior  resistance  for  particle  transport,  smaller  resistance  results 
in  higher  transport  rates,  which  indicates  that  the  surface  resistance  is 
dominant  in  the  process. 

The  variations  of  Re,  Sc,  and  Tr  on  the  steady  state  deposition  rates 
are  investigated  with  the  BLASIUS  program.  Figure  13  examines  the  Re 
number  variation  for  several  Tr  numbers  at  constant  Sc  of  10^.  On  a  log- 
log  scale  the  curves  are  all  fairly  straight.  They  more  or  less  fit  the 

relationship  of  Nu  «  Re^^.  The  role  of  Sc  number  is  examined  in 

s  .s . 

Figure  14.  The  straight  line  portion  of  the  curve  for  Tr=10^  is 

represented  by  Nu  <*  Sc  '  quite  well.  Also  Nu  increases  quickly 

s  •  s  •  s  •  s • 

.  —a  , 

when  Sc  is  larger  than  10  for  Tr=10  .  Figure  15  shows  the  effects  of  Tr 

number.  Based  on  the  variations  of  the  Nu  as  a  function  of  Tr, 

s .  s .  ’ 

Nu  saturates  when  Tr  is  smaller  than  10~^.  This  corresponds  to  the 

case  where  the  wall  is  a  perfect  sink  and  there  is  no  surface  resistance  to 
particle  transport  and  the  boundary  conditions  (52)  can  be  replaced  by 
(69).  As  Tr  gets  larger,  Nu  depends  more  heavily  on  its  variation. 

S  •  S  • 

CONCLUSION 

The  transport  and  deposition  of  fine  particles  in  a  boundary  layer  of 
impulsive  flow  over  a  flat  plate  have  been  examined  by  numerical  methods. 
The  transient  impulsive  flow  over  a  flat  plate  is  solved  by  a  vortex  sheet 
method  and  the  particle  transport  equation  is  numerically  integrated  by  the 
Strongly  Implicit  Procedure.  The  transient  solutions  are  preliminary 
because  refinement  of  the  computer  program  is  needed.  Some  trends  and 
interesting  features  may  be  drawn  from  the  results. 


The  transport  and  deposition  of  the  fine  particles  (d^<10u)  are 
;overned  by  Re,  Sc,  and  Tr .  For  smaller  values  of  Tr,  the  dimensionless 


deposition  rates.  Nu  is  proportional  to  Re*^  and  Sc  ^  ^  , 

s  •  s  • 


The  time  to 


reach  steady  state  seems  to  be  proportional  to  Sc,  inversely  proportional  to 
Tr  and  independent  of  Re. 

The  combination  of  a  vortex  sheet  method  for  the  flow  calculations  and 
a  finite  difference  method  for  the  transport  equation  should  be  examined 
more  closely.  The  vortex  sheet  method  is  designed  for  large  Reynolds  number 
flow  and  it  is  grid  free;  therefore,  the  solutions  can  contain  fluctuations 
due  to  statistical  errors.  Those  errors  may  be  reduced  by  averaging  the 
instantaneous  solutions.  By  applying  the  instantaneous  velocity  profiles 
directly  into  the  finite  difference  transport  equation,  the  statistical 
errors  will  amplify  the  truncation  errors.  It  is  thought  that  the  flow 
solutions  by  the  vortex  sheet  method  should  be  averaged  before  coupling  them 
with  the  transport  equation. 


RECOMMENDATIONS 

1.  While  the  instantaneous  flow  field  predicted  by  the  vortex  sheet 
method  represents  more  or  less  the  actual  nature  of  a  high  Reynolds  number 
flow,  a  Lagrange  type  of  approach  is  more  appropriate  for  investigating  the 
transport  of  larger  particles.  The  trajectory  of  each  particle  is  followed 
by  solving  the  equation  of  motion  for  this  particle.  A  large  number  of 
sample  particles  should  be  introduced,  then  the  overall  effects  are 
obtained  by  performing  statistical  analysis. 

2.  The  refinement  of  the  PARTRAN  computer  code  should  be  carried  out 
with  an  effort  first  to  investigate  how  to  reduce  the  statistical  errors  of 
the  random  vortex  sheet  method.  Next,  investigate  the  method  of  reducing 
the.  amplification  of  the  statistical  method  on  the  truncation  errors. 


SUMMARY 


In  many  two-phase  applications,  particle  deposition  in  a  boundary 
layer  is  strongly  dependent  upon  the  flow  field  and  the  particle  behavior 
near  the  wall  surface.  The  boundary  layer  is  a  shear  flow  with  different 
streamwise  and  vertical  characteristics  which  affect  the  particle  transport 
and  deposition.  In  transient  condition,  the  developing  impulsive  boundary 
layer  and  the  reverse  flow  will  influence  both  the  particle  transport  and 
deposition. 

The  present  study  is  concerned  with  the  transport  characteristics  of 
fine  particles  (<10p)  in  an  impulsive  boundary  layer.  A  random  vortex 
sheet  method  is  used  for  the  flow  generation  and  the  transport  equation  is 
solved  by  the  Strongly  Implicit  Procedure  of  finite  difference  scheme.  The 
particle  depositions  are  given  in  terms  of  the  Reynolds  number,  the  Schmidt 
number  and  the  transport  number  Tr  =  (2D/v+pL)  .  For  small  values  of  Tr, 
the  dimensionless  deposition  rate,  Nu  is  proportional  to  Re  ^  ^  and  to 

1  /O  * 

Scl,J.  The  time  to  reach  steady  state  seems  to  be  proportional  to  Sc, 
inversely  proportional  to  Tr  and  independent  on  Re.  In  general,  the 
deposition  rates  increase  with  increasing  Re  and  Sc  but  with  decreasing  Tr. 

The  combination  of  the  vortex  sheet  method  and  the  finite  difference 
method  at  the  current  form  described  in  the  report  does  not  predict  the 
steady  state  results  for  some  intermediate  Sc.  A  refinement  for  the 
program  is  needed  to  reduce  the  amplification  of  the  statistical  errors  due 
to  the  random  vortex  method  on  the  truncation  errors  of  the  finite 
difference  method  in  the  transport  equation. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the 
reports  it  publishes.  Your  comment s/answers  to  the  items/questions  below  will 
aid  us  in  our  efforts. 

1.  BRL  Report  Number _ Date  of  Report _ 

2.  Date  Report  Received _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or 

other  area  of  interest  for  which  the  report  will  be  used.) _ 


4.  How  specifically,  is  the  report  being  used?  (Information  source,  design 
data,  procedure,  source  of  ideas,  etc.) _ 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far 
as  man-hours  or  dollars  saved,  operating  costs  avoided  or  efficiencies  achieved, 
etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future 
reports?  (Indicate  changes  to  organization ,  technical  content,  format,  etc.) 


Name 


CURRENT 

ADDRESS 


Organization 


Address 


City,  State,  Zip 


7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the 
New  or  Correct  Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD 

ADDRESS 


Organization 
Address  " 


City,  State,  Zip 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


